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ORIGINAL ARTICLE 

Physiologically Based Pharmacokinetic Modeling 
Framework for Quantitative Prediction of an Herb-Drug 
Interaction 

SJ Brantley 1 , BT Gufford 2 , R Dua 1 , DJ Fediuk 1 ,TN Graf 3 ,YV Scarlett 4 , KS Frederick 5 , MB Fisher 6 , NH Oberlies 3 and MF Paine 2 

Herb-drug interaction predictions remain challenging. Physiologically based pharmacokinetic (PBPK) modeling was used to 
improve prediction accuracy of potential herb-drug interactions using the semipurified milk thistle preparation, silibinin, as an 
exemplar herbal product. Interactions between silibinin constituents and the probe substrates warfarin (CYP2C9) and midazolam 
(CYP3A) were simulated. A low silibinin dose (160mg/day x 14 days) was predicted to increase midazolam area under the curve 
(AUC) by 1 %, which was corroborated with external data; a higher dose (1 ,650 mg/day x 7 days) was predicted to increase midazolam 
and (S)-warfarin AUC by 5% and 4%, respectively. A proof-of-concept clinical study confirmed minimal interaction between high- 
dose silibinin and both midazolam and (S)-warfarin (9 and 13% increase in AUC, respectively). Unexpectedly, (f?)-warfarin AUC 
decreased (by 15%), but this is unlikely to be clinically important. Application of this PBPK modeling framework to other herb- 
drug interactions could facilitate development of guidelines for quantitative prediction of clinically relevant interactions. 
CPT Pharmacometrics Syst. Pharmacol. (2014) 3, e107; doi: 1 0.1 038/psp.201 3.69; published online 26 March 2014 



Herbal products represent an ever-increasing component of 
Western pharmacotherapy due in part to the perception that 
"natural" equates with safety. Current regulatory oversight 
of herbal products in many Western countries, including the 
United States, does not include evaluation of drug interac- 
tion liability prior to marketing. Consequently, systematic 
approaches for quantitative prediction of both the magnitude 
and likelihood of herb-drug interactions are nonexistent. 

Drug interaction liability assessment of herbal products is 
more challenging than for conventional drugs because unlike 
most drug products, herbal products typically are mixtures of 
bioactive constituents that vary substantially between prepa- 
rations. 1-3 Compounding this complexity is the often scant 
knowledge of specific causative constituents or systemic 
exposure of such constituents. Due to incomplete absorp- 
tion and extensive presystemic clearance, herbal product 
constituents may reach sufficient concentrations in the intes- 
tine and liver to inhibit only first-pass extraction of sensi- 
tive substrates. 4 Consequently, traditional (static) prediction 
approaches frequently do not translate to the clinical setting 
for herb-drug interactions. 

Recent drug-drug interaction guidelines suggest dynamic 
modeling and simulation approaches to predict complex 
interactions. 5 6 Extension of this approach to herb-drug inter- 
actions is a logical step to facilitate prospective evaluation of 
these interactions. As with drug-drug interactions, 78 physi- 
ologically based pharmacokinetic (PBPK) modeling may be 
used to improve in vitro to in vivo extrapolation of herb-drug 
interactions. Well-characterized herbal products are needed 
to develop a quantitative framework. 



Milk thistle preparations are top-ten selling herbal products 
in the United States. 9 The crude extract, silymarin, contains 
at least seven flavonolignans and one flavonoid. 10 The semi- 
purified extract, silibinin, contains roughly a 1:1 mixture of 
the flavonolignans silybin A and silybin B and represents 
an exemplar herbal product for initial model development. 
First, silybin A and silybin B have been purified in quanti- 
ties sufficient to recover requisite in vitro parameters. 1112 
Second, in vitro studies have demonstrated both reversible 
and mechanism-based inhibition of the key drug metaboliz- 
ing enzymes CYP2C9 13-15 and CYP3A4. 131416 Third, in vitro 
to in vivo extrapolation has been inconsistent. 17-19 Based on 
these observations, the objective of this work was to advance 
the mechanistic understanding of this herb-drug interaction 
using PBPK modeling and simulation with warfarin and mid- 
azolam as probe substrates. The models were evaluated 
through a proof-of-concept clinical study in healthy volun- 
teers. Results could help develop guidelines for prospective 
evaluation of herb-drug interaction liability. 

RESULTS 

Modeling and simulation 

PBPK model generation and initial evaluation. Simulated 
probe substrate concentrations closely approximated previ- 
ously published concentration-time profiles for both warfa- 
rin 20 and midazolam 21 under baseline conditions (data not 
shown). Model-predicted primary endpoints (area under the 
curve (AUC) and C max for (S)-warfarin and midazolam) were 
within the prespecified criterion (30%) for acceptable model 
performance (Table 1). 
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Table 1 Comparison of previously published and model-predicted 
pharmacokinetic outcomes 



Outcome 


Prpviou^lv 

published 3 


Model 
predicted" 


(%) 


(H)-Warfarin (5mg) c 








f 1/2 (hour) 


42(18) 


29 


69 


fnax (hour) (median 
(range)) 


2.0 (0.5-12) 


1.6 


80 


C max (umol/l) 


1 .7 (22) 


2.1 


124 


AUC 0 _ inf (umol/l-hour) 


93 (21) 


91 


99 


Cl/F (l/hour) 


0.18 (21) 


0.18 


100 


(S)-Warfarin (5mg) c 








t V2 (hour) 


32 (26) 


22 


69 


^max (hour) (median 
(range)) 


2 (0.5-4) 


1.5 


75 


C max (umol/l) 


2.0 (29) 


2.1 


105 d 


AUC 0 _ inf (umol/l-hour) 


65 (30) 


70 


108 d 


Cl/F (l/hour) 


0.25 (31) 


0.23 


92 


Midazolam (5mg) e 








t V2 (hour) 


2.9(41) 


3.5 


121 


f max (hour) (median 
(range)) 


0.5 (0.25-1.5) 


0.6 


120 


C max (nmol/l) 


88 (44) 


70 


80 d 


AUC 0 _ inf (nmol/l-hour) 


220 (33) 


210 


95 d 


Cl/F (l/hour) 


71 (33) 


72 


101 


Midazolam (8mg) f 








t V2 (hour) 


4.2 (29) 


3.5 


83 


f max (hour) (mean (SD)) 


0.47 (51) 


0.6 


128 


C max (umol/l) 


110 (49) 


110 


100 


AUC 0 _ inf (nmol/l-hour) 


300 (44) 


340 


113 


Cl/F (l/hour) 


95 (35) 


72 


76 


Silybin A (92.8mg)9 h 








t V2 (hour) 


1.6 


1.4 


— 


^max (hour) (median 
(range)) 


1 5 (1-2) 


1.3 




C max (umol/l) 


0.84 (89) 


0.27 


— 


AUC 0 _ 8 (umol/l-hour) 


1.3 


1.1 


— 


Cl/F (l/hour) 


150 


170 




Silybin B (128mg)9- h 








t V2 (hour) 


1.1 


1.4 




f max (hour) (median 
(range)) 


1 .5 (0.5-2) 


1.1 




C ma >mol/l) 


0.27(120) 


0.16 




AUC 0 _ 8 (umol/l-hour) 


0.28 


0.63 




Cl/F (l/hour) 


950 


410 





AUC 0 _ inf , area under the concentration-time curve from time zero to infinity; 
AUC 0 _ 8 , AUC from 0-8 hours; C max , maximal concentration; Cl/F, apparent oral 
clearance; f 1/2 , terminal half-life; f max , time to maximal concentration. 
a Geometric or arithmetic means and coefficients of variation (%) unless 
indicated otherwise. b Point estimates. ce_g Previously published outcomes from 
refs. 20,21 ,18,26, respectively. d Model predictions were considered accurate 
if the primary outcomes ((S)-warfarin and midazolam AUC and C max )) were 
within 30% of previously published outcomes. h Due to the sparse nature of 
the previously published data, a modified Bailer method (available in Phoenix 
WinNonlin) was used to recover t V2 , AUC 0 _ 8 , and Cl/F; accuracy was not 
calculated based on the sparse data and the 30% criterion being applicable 
only to the victim drugs. 




1 -I . 1 1 1 . » 

0 2 4 6 
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Figure 1 Mean concentration-time profile (0-6 hours) of midazolam 
in 19 healthy volunteers following an 8mg oral midazolam dose 
given alone (open symbols) or following a 1 4-day treatment with milk 
thistle product (solid symbols). 18 Lines denote physiologically based 
pharmacokinetic model simulations of the midazolam concentration- 
time profile when given alone (black) or with milk thistle (green). The 
dotted green line denotes incorporation of reversible inhibition of 
CYP3A, whereas the dashed green line denotes incorporation of 
mechanism-based inhibition of CYP3A. Symbols and error bars 
denote observed means and SDs, respectively, and were obtained 
from ref. 1 8. 

Prediction of silibinin-drug interaction magnitude. Simula- 
tions of a previously reported milk thistle-midazolam inter- 
action, assuming reversible CYP3A inhibition solely due to 
silybin A and silybin B, demonstrated negligible changes in 
the midazolam concentration-time profile (Figure 1). The 
milk thistle product tested, silymarin, contained 100mg 
of silybin A and 180mg of silybin B and was administered 
daily for 14 days. 18 Simulations assuming mechanism-based 
CYP3A inhibition predicted a 30% and 60% increase in mid- 
azolam C max and AUC, respectively; increases of 6% and 3% 
in midazolam C max and AUC, respectively, were reported 18 
(Figure 1). 

Simulations of the silibinin-warfarin interaction with a 
higher dose of silibinin (1 ,650mg/day, or 720 mg silybin A 
plus 930 mg silybin B/day; see below), assuming reversible 
CYP2C9 inhibition only, predicted negligible changes (<5%) 
in all pharmacokinetic outcomes (Figure 2a; Table 2). Sim- 
ulations of the high-dose silibinin-midazolam interaction 
assuming reversible CYP3A inhibition predicted no change 
in midazolam t V2 and <5% increase in both C max and AUC 
(Figure 2b; Table 2). Simulations assuming mechanism- 
based CYP3A inhibition predicted a 2-, 5-, and 1.5-fold 
increase in C max , AUC, and t V2 , respectively (Table 2). 

Proof-of-concept clinical evaluation 

Silibinin content in test product. A single lot (#304090) of 
Siliphos capsules, labeled to contain 60 mg silibinin, was 
selected. The capsules were overfilled consistently, contain- 
ing 69.1±4.28mg silibinin represented as 30.3±1.88mg 
silybin A and 38.9±2.39mg silybin B. The capsules also 
contained minor amounts of the regioisomers isosilybin A 
(1.55±0.09mg) and isosilybin B (0.94±0.06mg). 

Study subjects. All enrolled subjects completed the study 
(Supplementary Table S2). The study drugs and silibinin 
generally were well tolerated. One subject experienced mild 
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Figure 2 Geometric mean concentration-time profile of (a) warfarin, 
(b) midazolam, and (c) silibinin in 12 healthy volunteers following a 
10mg oral dose of warfarin or 5mg oral dose of midazolam given 
alone (open symbols) or following a 7-day treatment with silibinin 
(solid symbols). Lines in a and b denote physiologically based 
pharmacokinetic (PBPK) model simulations of the concentration- 
time profiles when the probe substrates were given alone (black) or 
with silibinin (green). Blue and orange lines in c denote PBPK model 
simulations of the concentration-time profiles of silybin A and silybin 
B, respectively. Symbols and error bars denote observed geometric 
means and upper limits of the 90% confidence interval, respectively. 

gastrointestinal upset following the first dose of silibinin. The 
effect was deemed likely to be drug related by the study phy- 
sician but was transient and did not limit the subject's contin- 
ued participation. During Clinical and Translational Research 
Center (CTRC) visits, four subjects (two in both phases) 
reported mild headaches attributed to caffeine withdrawal. 
No international normalized ratio elevations from baseline 
were observed following warfarin administration. 

CYP2C9 genotyping. All enrolled subjects consented to 
genotyping. Ten subjects were homozygous for the reference 
CYP2C9*1 allele. Two subjects carried one copy of the refer- 
ence allele and either the CYP2C9*2 or CYP2C9*3 allele. 

Effects of high-dose silibinin on warfarin and midazolam 
pharmacokinetics. The effects of high-dose silibinin 



(1,650mg/day) were compared to baseline oral pharma- 
cokinetics of warfarin and midazolam. Due to the reported 
mechanism-based inhibition of CYP3A in w'fro, 1416 silib- 
inin was administered three times daily for 6 days prior to 
administration of the probe substrates. Silibinin constituents 
were not expected to accumulate during the administration 
period due to short reported half-lives (<4 hours). 22 One 
subject demonstrated poor goodness-of-fit statistics for the 
t V2 of both warfarin enantiomers in both phases (FP < 0.85). 
Accordingly, data from this subject were excluded from 
analysis of AUC 0 48 and t V2 . 

Warfarin enantiomers were absorbed rapidly during both 
study phases, with median f max occurring at 1.25 and 1.5 
hours for (R)- and (S)-warfarin, respectively (Figure 2a). 
Coadministration with silibinin did not alter median (S)- 
warfarin f max but delayed median (ft)-warfarin f max by 15 
minutes. Relative to control (baseline), silibinin decreased 
(H)-warfarin geometric mean C max by 17% (Figure 3a; 
Table 2) and AUC 0 _ 48 by 15% (Figure 3b; Table 2). Silib- 
inin decreased geometric mean (S)-warfarin C max by 2% 
(Figure 3c; Table 2). Geometric mean AUC 0 48 of (S)-war- 
farin increased by 13% (Figure 3d; Table 2), with three 
subjects lying outside the predefined no effect range (0.75- 
1 .33). The 90% confidence intervals for the (S)-warfarin pri- 
mary endpoints (C max and AUC) lay within the predefined no 
effect range (Table 2). 

The rapid absorption of midazolam was unaltered by coad- 
ministration with silibinin, with median f max occurring at 0.5 
hours (Figure 2b). Relative to control, silibinin increased mid- 
azolam geometric mean C max by 20% (Figure 3e; Table 2) 
and AUC 0 inf by 9% (Figure 3f; Table 2). Except for one sub- 
ject (2.3-fold increase), treatment/control ratios of AUC 0 inf lay 
within the predefined no effect range (Figure 3f). The 90% 
confidence interval for midazolam treatment/control ratio of 
C extended above, whereas that of AUC n . , lay within, the 

max 0-inf J 7 

predefined no effect range (Table 2). 

The sampling strategy was not optimized for recovery of 
silybin A and silybin B pharmacokinetic outcomes; as such, 
these outcomes were interpreted for qualitative rather than 
quantitative purposes. The median f max of silybin A and sily- 
bin B following the initial administration of silibinin (3 and 3.5 
hours, respectively) nearly coincided with the second admin- 
istration of silibinin (Figure 2c). Geometric mean C max for 
silybin A was more than double that for silybin B (Table 2). 
Geometric mean t V2 of both silybin A and silybin B was ~5 
hours (Table 2). 



DISCUSSION 

Although herbal product usage continues to increase, cur- 
rent regulatory guidelines in several Western countries do 
not request premarket evaluation of herb-drug interaction 
liability. Investigations into such liabilities are fraught with 
inconsistent results due to the lack of a standard system 
for evaluation, high compositional variation between herbal 
products, and uncertainty about causative constituents. 
Unlike conventional drug products, the relative compo- 
sition of herbal products may vary substantially depend- 
ing on weather conditions, product collection and storage 
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Table 2 Comparison of proof-of-concept clinical study outcomes to physiologically based pharmacokinetic model predictions 



Observed 



Predicted (reversible inhibition) 



Predicted 
(mechanism-based inhibition) 



Outcome 


Geometric mean (CV %) 


Treatment/ 
control ratio (90% CI) 


Geometric mean 


Treatment/ 
control ratio 


Geometric mean 


- Treatment/ 


Control 


Treatment 


Control 


Treatment 


Control 


Treatment 


control ratio 


(H)-Warfarin 




















C max (umol/l) 


1 .92 (30) 


1 .60 (29) 


0.83 (0.77-0.90) 


2.08 


2.08 


1.00 











AUC 0 _ 48 (umol/l-hour) a 


55.0 (24) 


47.0 (23) 


0.85 (0.81-0.90) 


61.6 


61.5 


1.00 











t V2 (hour) 3 


52.0 (28) 


61.2 (27) 


1.14(0.96-1.36) 


29.0 


29.0 


1.00 











(S)-Warfarin 




















C max (umol/l) 


2.01 (32) 


1 .97 (27) 


0.98 (0.91-1 .05) b 


2.06 


2.08 


1.01 











AUC 0 _ 48 (umol/l-hour) a 


37.4 (41) 


42.3 (34) 


1.13 (1.01-1.26) b 


53.8 


56.0 


1.04 











f 1/2 (hour) 3 


29.6 (25) 


30.3 (20) 


0.97 (0.84-1.12) 


22.2 


22.4 


1.01 





, 





Midazolam 




















C max (nmol/l) 


74.2 (43) 


89.5 (39) 


1.20 (0.96-1 .51 ) b 


70 


73 


1.04 


70 


150 


2.11 


AUC 0 _ inf (nmol/l-hour) 


198 (42) 


216 (36) 


1.09 (0.93-1 .25) b 


210 


220 


1.05 


210 


1,070 


5.05 


f 1/2 (hour) 


5.17(36) 


4.90 (48) 


0.95 (0.82-1.10) 


3.55 


3.54 


1.00 


3.55 


5.30 


1.49 


Silybin A 




















C max (umol/l) 




0.97 (91) 






0.76 










f 1/2 (hour) 




5.1 (34) 






1.4 










Silybin B 




















C max (umol/l) 




0.40 (110) 






0.43 










f 1/2 (hour) 




5.1 (56) 






1.4 











half-life. 

a Evaluable for 1 1 subjects; all other outcomes were evaluable for 12 subjects. b The predefined no effect range was 0.75-1 .33 for the primary endpoints 
(S)-warfarin and midazolam treatment/control ratio for C max and AUC). 



methods, and processing procedures. 4 Accurate predic- 
tions of herb-drug interaction liability require not only iden- 
tification and quantification of causative constituents, but 
also measures of exposure in organs with metabolic capa- 
bility. Silibinin was selected as an exemplar herbal prod- 
uct due to a well-characterized composition, availability of 
inhibitory kinetic parameters from individual constituents, 
and disparate impact of milk thistle products on victim drug 
pharmacokinetics in previous clinical studies. 17-19 A PBPK 
modeling and simulation approach was used to address 
the challenges inherent to investigation of herb-drug inter- 
action liability. 

Warfarin is a widely used oral anticoagulant with a narrow 
therapeutic window. Warfarin is associated with a notoriously 
complicated pharmacotherapy due in part to myriad drugs 
and herbal products that alter the metabolism or anticoagu- 
lant activity of warfarin. As the clearance of the more phar- 
macologically active (S)-enantiomer is mediated primarily 
by CYP2C9, inhibition of this enzyme can lead to increased 
risk of bleeding. Silymarin was shown previously to increase 
systemic exposure to the CYP2C9/3A substrate losartan, 19 
prompting evaluation of the interaction potential between milk 
thistle and warfarin. Of the milk thistle constituents whose 
CYP2C9 interaction liability has been evaluated in vitro, sily- 
bin A and silybin B were the most potent. 15 These observa- 
tions led to the selection of silibinin, which consists primarily 
of these two constituents, for clinical evaluation. 



Relative to control, silibinin unexpectedly decreased 
both the geometric mean C max and AUC 0 last of (ft)-warfarin. 
Clinical manifestation of the previously reported CYP1A2 
induction by a milk thistle extract 23 is consistent with this 
decrease in exposure. In contrast to the doubling of losar- 
tan exposure following administration of silymarin, high- 
dose silibinin did not increase geometric mean (S)-warfarin 
exposure to a clinically relevant extent. However, increases 
above 33% were observed in three subjects, indicating 
that the CYP2C9 interaction potential of silibinin cannot 
be disregarded completely. Consistent with the expected 
decrease in (S)-warfarin clearance, the subject carrying 
the reduced function CYP2C9*3 allele demonstrated pro- 
longed warfarin exposure, which was not captured in the 
48-hour sampling window. As such, the AUC and t V2 for this 
subject were excluded from analysis. 

Modeling and simulation of the silibinin-warfarin interac- 
tion demonstrated that the rapid clearance of the silibinin 
constituents precluded marked inhibition of warfarin clear- 
ance. Sensitivity analysis of this interaction potential dem- 
onstrated that 10-fold increases in silybin A or silybin B 
inhibition potency (reversible K) would lead to roughly 15% 
increases in (S)-warfarin AUC (Supplementary Figure S1). 
Extensive intestinal and hepatic conjugation of silybin A and 
silybin B followed by rapid elimination likely would limit the 
interaction potential to first-pass clearance of sensitive sub- 
strates. Warfarin is not sensitive to first-pass elimination and 
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is cleared only upon subsequent passages through the liver, 
at which time any reversible inhibition of CYP2C9 by silybin 
A and silybin B would be abated. In contrast, losartan has a 
low bioavailability (33%) that is attributed, in part, to first-pass 
elimination. 24 This observation, coupled with the differences 
in study population and herbal product tested, could explain 
the difference between the reported interaction with losar- 
tan 19 and the lack of interaction with warfarin in the present 
study. Collectively, these observations suggest examination 
of other CYP2C9 substrates sensitive to first-pass elimina- 
tion, such as fluvastatin, to understand fully the milk thistle- 
CYP2C9 interaction potential. 

Midazolam is a gold standard CYP3A probe substrate 
metabolized extensively by intestinal and hepatic enzymes. 
Inhibition of CYP3A at either site can increase systemic 
exposure to midazolam; inhibition of hepatic CYP3A also 
can increase t V2 . Milk thistle constituents, including silybin A 
and silybin B, have been shown to be reversible and mecha- 
nism-based inhibitors of CYP3A activity in both human liver 
microsomes and expressed enzyme systems. 13 ' 1416 Previ- 
ous clinical interaction studies with midazolam 1718 have 
demonstrated limited interaction liability with the milk thistle 
product silymarin, albeit the doses administered were not 
sufficient to determine the difference between reversible 
and mechanism-based inhibition of CYP3A (Figure 1 ). The 
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Figure 3 Effects of silibinin (1 ,650 mg/day for 7 days) on (a,c,e) C max 
and (b,d,f) AUC of (a,b) (fl)-warfarin, (c,d) (S)-warfarin, and (ej) 
midazolam in 12 healthy volunteers following oral administration of 
warfarin (10mg) and midazolam (5mg). Open symbols connected 
by solid lines denote individual values. Solid symbols connected by 
dashed lines denote geometric means. 



"supratherapeutic" silibinin dose in the current study was 
selected to provide a large range between the predicted 
interaction based on reversible and mechanism-based inhi- 
bition of CYP3A and to maximize the ability to observe a 
clinical interaction. The lack of an interaction observed in all 
but one subject indicated that the CYP3A interaction liabil- 
ity for silibinin is low and is more consistent with revers- 
ible than mechanism-based inhibition (assuming inhibition 
indeed occurred). The current work represents another 
example of a potential mechanism-based inhibitor identi- 
fied in vitro™ that does not manifest clinically. 

Modeling and simulation of the silibinin-midazolam inter- 
action indicated that the low interaction potential is due, 
in part, to the lower inhibition potency of the silibinin con- 
stituents toward CYP3A compared to CYP2C9 (Table 3). 
Ten-fold increases in inhibition potency of silybin A and 
silybin B toward CYP3A activity increased midazolam 
exposure by roughly 25% (Supplementary Figure S1 and 
Supplementary Materials and Methods). These observa- 
tions indicated that at the predicted exposures, the constit- 
uents would need to be 10-fold more potent to demonstrate 

Table 3 Physiologically based pharmacokinetic model input parameters 



Victim compound 



Perpetrator 
compound 



Parameter 



(/?)- (S)- Silybin Silybin 

Warfarin Warfarin Midazolam A B 



Physicochemical/binding 



Molecular weight 


308.33 


308.33 


325.78 


482.44 


482.44 


Fraction absorbed 


1.0 a 


1.0 a 


1.0 a 


0.77 b 


0.77 b 


k a (hour 1 ) 


3.0 C 


3.0 C 


1.1 7 d 


0.50 b 


0.50 b 


Blood/plasma ratio 


1.0 e 


1.0 e 


0.80 f 


0.58 b 


0.58 b 


Unbound fraction in 
plasma 


0.0069 


0.0069 


0.02 f 


0.04 h 


0.04 h 


Metabolism 












Intestinal K m (umol/l) 






3.7 1 


22 c 


8.5 C 


Intestinal V max (umol/ 
hour) 






1,100' 


2,700 c 


2,600 c 


Hepatic K m (umol/l) 




6.5 e 


6.0 1 


54 c 


57 c 


Hepatic V max (umol/hour) 




260 c 


18,000' 


2,300 c 


2,700 c 


Hepatic Cl int (l/hour) 


30.4 C 










Inhibition 












CYP2C9 K. (umol/l) 


6.5J 






10 k 


4.8 k 


CYP2C9 a 








5 k 


8 k 


CYP3A4 K. (umol/l) 








26.5' 


31 .5' 


CYP3A4 k inact (minute- 1 ) 








0.22' 


0.15' 


CYP3A4 K, (umol/l) 








100' 


89' 



See "Methods" for detailed information on model parameterization. 
Cl int , intrinsic clearance; k a , absorption rate constant; K., reversible inhibition 
constant; a, affinity change of the enzyme-substrate and enzyme-inhibitor 
complexes; K v concentration required to achieve half-maximal rate of enzyme 
inactivation (/c. nact ). 

a Assumed. Predicted based on physicochemical properties using ADMET 
Predictor (Simulations Plus). c Obtained by fitting the model to clinical data 
(ref. 26). d Ref. 35. e Ref. 34. f Ref. 30. sRef. 33. h Ref. 40. Extrapolated from in 
vitro data. j Ref. 34. k Ref. 15. 'Obtained from recombinant data in ref. 24. 



www.nature.com/psp 



PBPK Modeling Framework for Quantitative Prediction of an Herb-Drug Interaction 

Brantley eta/. 



any clinically relevant interaction with CYP3A. The large 
predicted increase in midazolam exposure incorporating 
mechanism-based inhibition further supported the hypoth- 
esis that products with limited systemic exposure (first 
posited with fruit juices) 25 need to be mechanism-based 
inhibitors of CYP enzymes to perpetrate clinically relevant 
interactions. 

One limitation to the current work is that silybin A and 
silybin B clearance parameters were recovered by fitting 
the model to data obtained from hepatitis C patients admin- 
istered a product (silymarin) that contained additional con- 
stituents not present in silibinin. 26 In vitro determination of 
silibinin clearance parameters would provide a true bottom- 
up modeling approach and reduce complexities inherent to 
pharmacokinetic data from patients with hepatic disease. 
Alternatively, disease-related parameters could be used to 
develop a hepatitis C virtual population before fitting the 
PBPK model with the observed pharmacokinetic data, facil- 
itating recovery of disease-independent silibinin clearance 
parameters. 27-29 

In summary, prospective evaluation of herb-drug inter- 
actions, consistent with that for drug-drug interactions, 
largely has been ignored due to substantial compositional 
variability inherent to herbal products, multiple inhibitory 
constituents, varying inhibition mechanisms, and relative 
lack of regulatory oversight. The PBPK interaction model 
developed in the current work incorporated in vitro inhibi- 
tion kinetic parameters and systemic exposure estimates 
of individual constituents for the exemplar herbal product, 
silibinin. Simulations of the silibinin-warfarin and silibinin— 
midazolam interactions accurately predicted minimal clini- 
cal interaction liability. This work demonstrated the utility 
and predictive power of PBPK modeling and simulation, 
which could be extended to investigate scenarios (e.g., 
wide dosing ranges, tissue exposure assessment, and 
herbal product composition variation) and patient popula- 
tions (e.g., pediatric, geriatric, and pregnant women) not 
amenable to clinical investigation. Refinement of the PBPK 
model by recovering disease-independent silibinin clear- 
ance parameters and incorporating alternate victim drugs, 
including losartan, will enhance confidence in model pre- 
dictions and generalizability. This framework represents an 
initial step to establishing a systematic approach that can 
be applied to other combinations of herbal products and 
conventional drugs under various clinical scenarios to iden- 
tify potential clinically significant herb-drug interactions, 
predict the extent of those interactions, and ultimately help 
guide pharmacotherapeutic decisions. 

METHODS 

PBPK model development The base model structure 
was adapted from the literature 30 (Figure 4), incorporat- 
ing physiologic parameters obtained from the International 
Commission on Radiological Protection. 31 Warfarin partition 
coefficients (K p s) 32 and binding parameters 33 were obtained 
from the literature (Table 3); absorption rate constants (k a s) 
and clearance parameters were obtained by fitting the PBPK 
model to previously reported plasma concentration-time 




Figure 4 Base physiologically based pharmacokinetic model 
structure. Model structure was modified from the literature. 30 Organ 
weights and blood flows were obtained from the International 
Commission on Radiological Protection. 31 Following oral 
administration, drug transfer from dosing compartment to intestine 
is driven by the oral absorption rate constant (k a ). Drug clearance 
(CI) is mediated by metabolic processes in the intestine and liver. 
The pancreas and spleen were combined into a hybrid "organ" 
designated as PSP. 

profiles. 20 The reversible inhibition constant (K) of (R)- 
warfarin toward CYP2C9 activity was obtained from the 
literature. 34 Midazolam K s and k were obtained from the 

P a 

literature 3035 ; intestinal and hepatic clearance parameters 
were extrapolated from in vitro data 16 as described 3637 
(Table 3). Silybin A and silybin B K p s were predicted from 
physicochemical properties 38 using GastroPlus (version 8.0; 
Simulations Plus, Lancaster, CA). Silibinin binding param- 
eters were obtained from the literature 39 ; clearance param- 
eters were generated by fitting the PBPK model to plasma 
concentration-time data from hepatitis C patients receiving 
silymarin 26 (Table 3). Silybin A and silybin B mechanism- 
based (Kj, /c nact ) and reversible inhibition kinetic parameters 
were obtained from the literature. 1516 Mechanism-based 
inhibition of CYP2C9 was not considered based on a previ- 
ous publication showing no IC 50 shift using (S)-warfarin as 
the probe substrate. 15 

PBPK interaction model simulations. PBPK models were 
developed for midazolam, (ft)-warfarin, (S)-warfarin, silybin 
A, and silybin B using Berkeley Madonna (version 8.3; Uni- 
versity of California at Berkeley, Berkeley, CA) with code 
compiled in MEGen 40 (version 0.5; UK Health & Safety 
Laboratory, Buxton, UK) (Supplementary Materials and 
Methods). The PBPK model for perpetrator (silybin A and 
silybin B) and victim (warfarin or midazolam) compounds 
were linked through the reversible or mechanism-based 
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inhibition of victim probe substrate. Initial simulations used 
doses of probe substrates and milk thistle products reported 
in previous studies. Simulations were considered accurate 
if the predicted primary pharmacokinetic outcomes (AUC 
and C max for (S)-warfarin and midazolam) were within 30% 
of observed outcomes. Following initial model evaluation, 
simulations were conducted with a higher dose of silibinin 
(1,650mg/day) to determine whether a clinically important 
interaction is possible. Pharmacokinetic outcomes from the 
simulated profiles were recovered via noncompartmental 
analysis using Phoenix WinNonlin (version 6.3; Pharsight, 
Cary, NC). 

Analysis of silibinin product. Siliphos capsules (n = 28) 
(Thorne Research, Dover, ID) were analyzed using a modi- 
fication of previously described methods 4142 to ensure purity 
and content. Briefly, the contents of each capsule were 
weighed and extracted twice with 2 ml acetone. The extract 
was vortex mixed and centrifuged (13,000gfor 2 minutes); 
the supernatant was transferred to a clean vial. Milk thistle 
constituents were quantified using an Acquity UPLC system 
with an HSS-T3 1.8 urn (2.1 x 100 mm) Acquity column and 
Empower 3 software (Waters, Milford, MA). Standards and 
Siliphos capsule extracts were analyzed using a gradient 
from 30:70 to 55:45 methanokwater (0.1% formic acid) over 
5.0 minutes at a flow rate of 0.6 ml/minute at 50 °C; peaks 
were detected at 288 nm. 

Proof -of -concept clinical study. Healthy volunteers (six men 
and six nonpregnant women) were enrolled in an open- 
label, fixed sequence crossover study conducted at the 
UNC CTRC. The study protocol was approved by the UNC 
Office of Human Research Ethics Biomedical Institutional 
Review Board and the CTRC Advisory Committee. Eligibility 
to participate was based on screening evaluation and inclu- 
sion/exclusion criteria (Supplementary Table S1). Written 
informed consent was obtained from each subject prior to 
enrollment. 

The first (control) phase consisted of administration of 
10mg warfarin (Coumadin; Bristol Meyers Squibb, Princ- 
eton, NJ), 10mg vitamin K (Mephyton; Aton Pharma, 
Lawrenceville, NJ), and 5mg midazolam syrup (Ranbaxy; 
Jacksonville, FL). A negative pregnancy result was required 
before drug administration to women of childbearing poten- 
tial. Vital signs (blood pressure, temperature, respiratory 
rate, pulse, and oxygen saturation) were obtained at base- 
line and every 15 minutes for the first 2 hours. All subjects 
underwent an international normalized ratio with prothrom- 
bin time. Blood (7 ml) was collected through an intravenous 
line before and from 0.25-12 hours following drug admin- 
istration. Subjects continued to fast until after the 4-hour 
blood collection, when meals and snacks, devoid of fruit 
juices and caffeinated beverages, were provided. Subjects 
returned to the CTRC 24 and 48 hours post-drug administra- 
tion for blood collection. Optimal study design simulations 43 
of previously reported clinical data 20 demonstrated that a 
0-48-hour collection was an accurate surrogate of total sys- 
temic exposure (AUC 0 inf ) for warfarin. Plasma was collected 
and stored at -80 °C pending analysis by high-performance 



liquid chromatography tandem mass spectrometry (HPLC/ 
MS-MS). 

Following at least a 14-day washout, subjects received 
480 mg silibinin (based on labeled content) to self-adminis- 
ter three times daily for 7 days. Each subject received his/ 
her silibinin in a blister pack and was asked to complete a 
pill diary documenting the time of administration. Subjects 
were contacted at least twice during the week of silibinin self- 
administration to monitor compliance and adverse events, 
which were graded using a validated Adverse Events Scale. 
Subjects returned to the CTRC on day 7 for concomitant 
administration of silibinin, warfarin, vitamin K, and mid- 
azolam. Plasma was collected and stored as described for 
the first phase. 

Analysis of plasma for warfarin enantiomers, midazolam, 
silybin A, and silybin B. Concentrations of all analytes were 
quantified using a Sciex (Framingham, MA) API4000 Qtrap 
HPLC-MS/MS triple quadrupole mass spectrometer fitted 
with a Turbo ionspray interface operated in the positive ion 
mode. Plasma was treated with acetonitrile (6 volumes) 
containing the internal standard, warfarin-d 5 (Toronto 
Research Chemicals; Toronto, Canada) or I'-hydroxymid- 
azolam-d 4 (Cerilliant, Round Rock, TX), and centrifuged 
(3,000g). The supernatant was injected into the HPLC-MS/ 
MS system. Warfarin enantiomers were separated on a 
Supelco Astec Chirobiotic V 15cm x 2.1 mm 5 micron chi- 
ral column (Sigma Aldrich; St Louis, MO) and eluted with 
an isocratic mixture consisting of 75% 5 mmol/l ammonium 
acetate containing 0.01% (v/v) formic acid and 25% aceto- 
nitrile (flow rate, 0.4 ml/minute). Midazolam was eluted with 
a binary gradient mixture consisting of 10 mmol/l ammo- 
nium formate containing 1% (v/v) isopropyl alcohol and 
0.1% (v/v) formic acid and methanol on a Varian Polaris 
C18-A 20cm x 2.0mm 5 micron column (Agilent, Santa 
Clara, CA) (flow rate, 0.65 ml/minute). Silybin A and silybin 
B were eluted with an isocratic mixture consisting of 44% 
water, 56% methanol, and 0.1% (v/v) formic acid on an Agi- 
lent Zorbax XDB C18 15cm x 3.0mm 3.5 micron column 
(Agilent) (flow rate, 0.7 ml/minute). Analyte concentrations 
were determined by interpolation from a linear standard 
curve with an assay dynamic range of 0.5-10,000 nmol/l 
(warfarin enantiomers) or 0.5-5,000 nmol/l (midazolam, 
silybin A, silybin B). Analytical methods were validated 
according to US Food and Drug Administration guide- 
lines. 44 Inter- and intraday variability for all analytes was 
less than 10%. 

Pharmacokinetic analysis. Pharmacokinetic outcomes were 
recovered by noncompartmental analysis using Phoenix 
WinNonlin. Concentrations below the limit of quantification 
were excluded. The terminal elimination rate constant (k z ) 
was estimated by linear regression of the terminal portion of 
the log-transformed concentration-time profile using at least 
three data points. The terminal half-life (t V2 ) was calculated as 
\n2/X . The maximum observed concentration (C ), time to 

z v max' ' 

reach C (t), and last measured concentration (G) were 

max v max' ' x last' 

obtained directly from the concentration-time profile. AUC 
from time zero to C last (AUC 0 last ) was determined using the 
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trapezoidal method with linear up/log down interpolation. The 
AUC from time zero to infinity (AUC 0 inf ) was calculated as the 
sum of AUC 0 last and the ratio of C last to X z . 

Genotyping for common CYP2C9 variants. CYP2C9*2 and 
*3 polymorphisms were determined using a previously pub- 
lished polymerase chain reaction restriction fragment length 
polymorphism assay. 45 

Statistical analysis. All statistical analyses were conducted 
using SAS (version 9.2; SAS Institute, Cary, NC).The sample 
size for the proof-of-concept study (n = 12 evaluable sub- 
jects) was calculated based on 80% power to detect a 25% 
change in the primary endpoints with a type I error of 0.05; the 
primary endpoints were the treatment/control ratios of log- 
transformed AUC 0 48 ((S)-warfarin) or AUC 0 inf (midazolam) 
and C max ((S)-warfarin and midazolam), and the predefined 
no effect range was 0.75-1 .33. 56 Intraindividual variability in 
midazolam and warfarin AUC and C were assumed to be 

max 

~20%. 46-48 Secondary outcomes, t V2 and f max , were evaluated 
using a paired two-tailed Student's ttest on log-transformed 
data or Wilcoxon signed-rank test as appropriate, with 90% 
confidence intervals and ranges reported for t V2 and f max , 
respectively. A P value <0.05 was considered statistically 
significant. 
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Study Highlights 

WHAT IS THE CURRENT KNOWLEDGE ON THE 
TOPIC? 

S Despite increasing recognition of herb-drug in- 
teractions in clinical practice, robust information 
about the causative ingredients and mecha- 
nisms underlying these interactions remains 
limited. Consequently, evidence-based recom- 
mendations about adding herbal products to ex- 
isting pharmacotherapeutic regimens virtually 
are nonexistent. 



WHAT QUESTION DID THIS STUDY ADDRESS? 

S This study addressed the utility of a PBPK 
modeling approach to predict the drug interac- 
tion liability of an herbal product. This approach 
was tested using the exemplar herbal product 
silibinin and the widely used cytochrome P450 
probe substrates warfarin and midazolam. 



WHAT THIS STUDY ADDS TO OUR KNOWLEDGE 

S A PBPK modeling approach accurately predict- 
ed the minimal interaction potential of chronic 
exposure to high-"dose" silibinin and two FDA- 
recommended probe substrates. Sensitivity 
analysis demonstrated that silibinin constituents 
are cleared too rapidly to influence the systemic 
metabolism of warfarin and that the inhibitory 
potency toward CYP3A is not sufficient for clini- 
cal interactions with midazolam. 



HOW THIS MIGHT CHANGE CLINICAL 
PHARMACOLOGY AND THERAPEUTICS 

S A PBPK modeling and simulation approach 
could facilitate prospective evaluation of herb- 
drug interactions, as well as evidence-based 
recommendations about adding herbal prod- 
ucts to conventional drug regimens. 
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